Mon. Not. R. Astron. Soc. 000, [T>?? (-) Printed 1 September 2010 (MN Wl&i style file v2.2) 



Corotating light cylinders and Alfven waves 



K.N. Gourgouliatos 1>2 * and D. Lynden-Bell 1 * 

1 Institute of Astronomy, University of Cambridge, Madingley Road CB3 OH A 

2 Department of Physics, Purdue University, 525 Northwestern Avenue, West Lafayette, IN 47907-2036, USA 



Accepted -. Received -; in original form - 



ABSTRACT 

Exact relativistic force free fields with cylindrical symmetry are explored. Such helds 
are generated in the interstellar gas via their connection to pulsar magnetospheres both 
inside and outside their light cylinders. The possibility of much enhanced interstellar 
fields wound on cylinders of Solar system dimensions is discussed but these are most 
likely unstable. 
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1 INTRODUCTION 



Despite the pioneering work of iGoldreich" fc Julian! (|l969l ) 
on pulsar mag netosph e res cr eated by a rotating aligned 
dipole and of IMichell j 19731 ) on relativistic winds both 
of which have been further explored in iMestell l|l999l ). 
there are few exact analytic solutions in relativistic mag- 
netohydrodynamics (MHD). Prendergast's (2005) solu- 
tion for a point explosion of force-free magnetic energy 
was the foundation of our recent work on that theme 
(IGourgouliatos fc Lvnden-Belll |2008l ; IGourgouliatos! 120091 ; 
IGourgouliatos fc VlahakisT 201C ) but recent advances have 



been mainl y computational. The earl y axisy mmetric static 
solution by IScharlemann fc Wagoner] (1 19731) has been fol- 



lowe d by accurate compu tations of Komissarovl (|2006t ) 
and IGoodwin et al.l (|2004l) which give solutions for the 
steady aligned pulsar while iKalapotharakos fc Contopoulosl 
(2009) considered both time-dependent instabilities at the 
edge of th e coro tating zone and the non-aligned case. 
ISpitkovskvl (|2006t ) studied the time dependent relativis- 
tic problem of both the aligned and the oblique ro- 
tator; in the case of the aligned rotator h e confi rmed 
the results of IContopoulos. Kazanas fc FendtJ fl999) who 
have studied the time-independent problem, whereas for 
the oblique rotator he presented the three-dimensional 
structure of the magnetosphere. Recent evidence that the 
pulses are gener ated near the light cylinder rein force ideas 
put forward by Ida Costa fc Kahnl (|l98l Il985l ) who con- 
sidered the orthogonally aligned cylindrical pulsar. Nu- 
merical solutions more rele vant to active galactic nu- 
clei h ave b een developed by iGammie. McKinnev fc Tothl 
J2003h and iMcKinnev fc Naravanl (120071). More recently 
iTchekhovskov. McKinnev fc Narava nl j2009l) have studied 
the problem of a spinning monopole, and they found that 
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the combined effect of rotation and the magnetic field pro- 
vide sufficient conversion of magnetic energy to kinetic with 
profound applications to GRBs but also to pulsar systems. 

Force-free fields have been used as an approximation for 
pulsar magnetospheres. Leaving a magnetic field configura- 
tion to relax, we take a final state of a force-free magnetic 
field. Taking into account that the pulsar magnetosphere 
is dominated by the electromagnetic fields and considering 
that the effects of inertia are weaker makes the force- free as- 
sumption a valid starting point. The non-relativistic force- 
free approximation requires that the electric current flows 
along the magnetic field lines but neglects any effect of the 
electric field on the charges. In relativistic MHD the elec- 
tric field is no more negligible as the electric field, in non- 
resistive MHD, is given by the cross product of the magnetic 
field with the velocity scaled to the speed of light. Thus the 
relativistic force-free approximation requires that the force 
of the magnetic field on the currents and that of the induced 
electric field on the charges add to zero. 

Here we explore the simplest exact relativistic force-free 
fields, those with cylindrical symmetry, and their coupling 
to rotating magnetospheres which corotate out to the light 
cylinder. While these solutions are very limited they allow 
us to answer questions such as: When one end of a cylinder 
of force-free field is rotated at relativistic speed, is its flux 
per unit area greater than or less than the static field that 
surrounds it in pressure equilibrium, i.e. if we start rotat- 
ing one end of a tube of flux, does the twist generated in 
it lead to its expansion or does its cross-section remain the 
same or even contract? What if anything happens as the 
rotation velocity of the tube approaches the speed of light? 
How rapidly does the rotation advance along the tube? For 
non-relat i vistic MHD this problem was treated in detail in 
ISturrockl (|l994h but inevitably that has little bearing on 
light cylinders. We are studying idealized MHD in which 
the effects of resistivity are absent, therefore field lines that 
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join the source of twist rotate while those that do not are left 
static. Thus a discontinuity is no surprise. In Section 2 we 
give the basic equations and solve them to find the general 
steady cylindrical force-free fields. Section 3 explores some 
simple special cases in greater detail, while Section 4 gives 
the time-dependent exact solution for rotationally generated 
Alfven waves traveling down a uniform field. Section 5 con- 
siders how a rotating star or a pulsar might generate such 
waves in the interstellar magnetic field. It also presents the 
conundrum of how a pulsar's field that extends beyond the 
light cylinder can connect to an anti-aligned interstellar field 
and suggests a possible resolution. 



known, E follows from ([3]), p from Q and j from ((5|. Thus 
the whole electromagnetic configuration is determined but 
to find the flow we must also specify the velocity Vb(R) 
along B. This is subject to the restriction that V 2 ^ 1 
which gives (RSl/c + VbB^/B) 2 + V£B 2 /B 2 sC 1 that is 
V£ + 2V B {RSl/c){B 4> /B) < 1 - R 2 Sl 2 /c 2 , which can be sat- 
isfied by a real Vb for example Vb = —2—^- provided 
1~("T^) 2 ("b l ) 2 ^0 which is always satisfied since from (J9J 
the expression on the left is < B 2 > / B 2 . Notice there is no 
restriction that RSl/c ^ 1 because the fluid can move along 
the lines of the force to counteract such a 'superluminal' ro- 
tation leaving its total V < 1. This is true however large 
RSl/c becomes! 



2 BASIC EQUATIONS 



We take the magnetic field and the velocity of the fluid to lie 
on cylinders and in cylindrical polars (R, 4>,z) so V • B = 
and 



B = [0,B^R),B z (R)}, 

v = cV = c[0,V^(R),V s (R)]. 

The fluid is perfectly conducting so 

E = V x B = [B V Z - fl„V*,0,0] 



(1) 
(2) 



(3) 



since Er is independent of <f> and z it follows that V x E = 
so B is time independent. The charge density, p, is given by 

V-E = 4irp=±JL\R(B 4 Vz-B z Vt)]. (4) 

Provided V is independent of time the current density, j, is 
given by 

4*3=VxB = -£ i B s $+±£ i (RB4)i. (5) 
The force-free condition is 

pE + j x B = . (6) 
Using the expressions for p and j above, this becomes 

^[R(V x B)f - \ ±Bl - ^(RB,) 2 = 0. (7) 

We notice that only the component of V perpendicular to B 
plays any role in equation (0 . It is then convenient to write 
V in the form V = c _1 fi(.R) x R + V B {R)B, where B = 
B/\B\ is the unit vector along B. We note that Vb is not 
the component of V along B which contains also (ft x R) ■ 
B/\cB\. Thus any of our velocity fields may be considered as 
a (differential) rotation with flow along B superposed. Since 
the the flow along B does not contribute to the (V x B) 2 
term in (0 that equation becomes 

Integrating the last term by parts and then dividing by R 2 
we find 



B\ =< B 2 > ~B 2 (1 - (c^SlR) 2 



(9) 



where < B 2 >= R~ 2 J Q B 2 dR 2 is the area average of B, 
up to R. Evidently if Q(R) and B Z (R) are specified such 
that the right hand side of ([9]) is non-negative, then ([9]) 
determines B^, (up to a sign). With B and V x B then 



3 SIMPLE SPECIAL CASES 

Looking at equation ((9} there are several special cases that 
are worth noticing because they illustrate situations of par- 
ticular simplicity. All the fields so far discussed are heli- 
cal so it is useful here to define the pitch k(R) of a heli- 
cal field line. The equation of such a field line is given by 
Rdt\>jB^ = dz/B z with R constant. The pitch is defined as 
k(R) = d<f>/dz so kR = B^/B z . 

(i) If B z is independent of R then < B 2 Z >= B z so from 
® B\ = {QR 2 /c 2 )B 2 . When SIR = c, B^ = ±B Z but as 
pointed out above there is no restriction that SIR < c so we 
can have situations with \B^,\ > \B Z \ at large R. There is 
however a tendency for such solutions to be unstable. In all 
cases the lines of force are given by R =const and = 

■g 2 -. The pitch of the resultant helix is k — d(j>/dz so kR = 
B^,/B z . 

(ii) Writing equation (|9]) in terms of the pitch k(R) we 

find 



B 2 z [1 + {k 2 -%)R 2 ] =<B 2 >, 



(10) 



A remarkable special case occurs when k 2 — S1 2 /c? is indepen- 
dent of R. This occurs when for example both the pitch and 
the angular velocity are both constant. Then we may mul- 
tiply equation (|10|) by R 2 and differentiate the result with 
respect to R 2 ; this gives us B z on the right which cancels 
against a term on the left. Dividing the resultant equation 
by R 2 yields 



dBl a SI 2 

r[l + (k 



)R 2 ]+2(k, 2 - Sl 2 /c 2 )B 2 



0. 



dR 2 ' C 
This integrates to give, on taking the square root 



\B*(R)\ = 



BJfi) 



1 + (k 2 -S1 2 /c 2 )R 2 



(11) 



(12) 



Notice that the field is uniform if k 2 = Sl 2 /c 2 . If now a 
rotating tube of force is surrounded by a static uniform field 
Bo at equilibrium with it at R — a 



Bl + B 2 + ATvaE R 



Bq , 



(13) 



where a is the surface charge on the rotating cylinder and Er 
is the electric field component inside the boundary R — a. 
By Gauss's theorem 4ira = —Er(cl) and from Q Er = 
{<r x SlaB*f. Hence \B z (a)\ = B /F 1/2 where F = 1 + (k 2 ~ 
Sl 2 /c 2 )a 2 . Thus when k — Sl/c the rotating uniform field 
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has the same \B Z \ as the non-rotating field outside. More 
generally the flux along the cylinder of radius a is 

J^2ivB z RdR = ttB z (0)(k 2 - n 2 /c 2 ) _1 lnF 
= vrB z (a)(K 2 - Q 2 /c 2 )- 1 F\nF 
= Tra 2 B [{F -l)~ 1 F 1 / 2 \nF]. 



(14) 



This function of F maximises at Tva 2 Bo when F = 1 that is 
when £7 = ck. In all other cases the equilibrium twisted field 
configuration carries less flux over the circle R < a than the 
straight bounding field Bo- Thus although some feel that 
the twisting of the field should cause it to contract under 
the influence of the hoop stresses so generated, nevertheless 
this is more than offset by the increase in field strength per 
unit longitudinal flux which is also caused by the twisting 
of the field. Notice however we have only proved this when 
k 2 — Q 2 /c 2 is independent of R within the twisted region. 
However for (k 2 — Q 2 /c 2 )a 2 between —0.5 and +1.0 the flux 
decrease is less than 2% so there is only a < 1% increase in 
the radius of the twisted rotating cylinder in such circum- 
stances. An untwisted B field spinning about the axis of a 
cylinder whose top speed is 1 /s/2 of the velocity of light only 
increases its radius by about about 1%! 

(iii) Now suppose that the rotation is uniform out to 
R = b but thereafter takes the Keplerian profile Q — 
Qo(b/R) 3 ^ 2 . We substitute this velocity profile in equa- 
tion © 

Bl=<B 2 z > -B 2 {l-n 2 R 2 /c 2 ) (15) 

Now the Alfven waves generated by such a rotation will 
travel at the speed light in a force- free field (see section 4) so 
the pitch generated is given by k, — d(j>/dz = —Q(R)dt / '(cdt) 
so B,p = B z Q(R)R/c so the first and last terms of equa- 
tion (I15p cancel so we are left with < B z >= B 2 so the 
z-component of the field is uniform despite its differential 
rotation i.e. we are back to (i). 



4 TIME DEPENDENT SOLUTIONS 

Since the solutions that leave the flux through the cylin- 
der unchanged have a uniform B z we now look for time 
dependent solutions of the form B — [0, B^(R, z, t), Bo] so 
V-B = 



v = cV = [0, Q(R, z, t)R, 0] + VbB , 
then, 

E = V x B = --B R, 

c 



Vx£= -B Rc~ l ^4> = --B<t,4>, 
az c 



1 • Q 

—E — — -BqR, 

c c/ 



(16) 

(17) 
(18) 
(19) 
(20) 
(21) 



c c z az 



}R + 



+ ilR^ = ^- 
From the force-free condition © we have 



(22) 



R 



d 



(bIr 2 



Bl d (SfR 4 

,2 



2R 2 dR 



n 



(23) 



The z and <f> components of this equation are satisfied pro- 
vided 



c l az 
but by (HHJ 



B R 



dz 



Bi, 



hence 

4« = 



d 2 n 

dz 2 



(24) 



(25) 



(26) 



which wave equation is satisfied by £1 = Q[R,(t — z/c)]. 
Then B<f, — —BqR£1/c and putting this back into (123[l the 
R component is satisfied identically, so we have an exact 
solution for any wave form f2(i — z/c). Of course we could 
superpose any £l(t + z/c) but to solve the problem of a per- 
fectly conducting disc at z — forced to rotate at a given 
time dependent rate Q(t) we need the wave that travels out 
from the origin. Then the magnetic field at z > 0, R < a 
rotates at Q(t — z/c) and that at z < 0, R < a rotates at 
Q(t + z/c), see Fig. 1. Thus we have Alfven waves gener- 
ated by a rotating disc. At R = a there is a surface density 
of charge 4tv<j = ^flaBo; fla < c, which accounts for the 
electric field outside R — a. The changing surface density 
of charge involves a sheet current along the cylinder given 
by integrating !)23|) across the discontinuity in B$ at R = a 
thus 47rJ z = B(,aQ./c and ^ aL + f- = as required by charge 
conservation. Note that the net charge of the rotating cylin- 
der is zero, as the volume charge density given by (Q and 
the surface charge density when integrated over the cylin- 
drical volume and the discontinuity surface respectively add 
up to zero, s o the electric field outside the cylinder is zero. 
iKomissarovl (|2002l 'l in the context of Blandford-Znajek mech- 
anism has investigated similar solutions. 



5 ROTATING MAGNETIC FIELDS IN THE 
INTERSTELLAR GAS 

So far our solutions have been exact but abstract. We now 
take a step nearer to the real world but the application is 
somewhat forced and the solution inexact. 

Consider a star with a dipolar magnetic field with axis 
along Oz. Let it be embedded in an interstellar magnetic 
field oriented in the same direction. Then the field is as 
shown in Fig. 2 where the separatrix has been emphasised. 
The closed field lines near the star will lie in meridional 
planes and rotate with the star. The outermost field lines 
outside any separatrix will not rotate at all just like the 
outer field lines in Fig. 1. The remaining field lines that join 
the star to the interstellar gas will carry the torque that 
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Figure 1. The magnetic field lines attached to a disc which is 
rotating rigidly. The twisted field is confined by a uniform exter- 
nal field. The heavy dashed line is the separatrix which carries a 
surface charge. When the tip speed of the disk is c this is also the 
light cylinder. 



generates the Alfven wave that propagates along the inter- 
stellar field to great distances from the star. In our force-free 
approximation these Alfven waves travel at the velocity of 
light but in reality the fields will probably be too weak well 
away from the star so there the rotational disturbance will 



travel at the Alfven speed given by (? A — 



where 



po is the gas density. In Fig. 2 we have assumed that the 
light cylinder lies beyond the asymptotic radius of the sep- 
aratrix. For a strong stellar dipole or a weaker interstellar 
field this will not be true. We shall consider that situation 
afterwards. We write the field in the form 



B = i(VP x V0 + /3V0) . 

Z7T 



(27) 



Then for a dipole M in a straight external field we find 



P = Tv{B r 2 +2Mr _1 )(l 



(28) 



The separatrix P = Ps — SitBq' 3 M 2 ^ 3 which for large z 
gives R s = a= V3(M/B ) 1/3 . Thus all the flux with P < 
Ps links the star to the interstellar magnetic field. At large 
z we have from Section 3 (i) and under equation (27) /3 = 
—2tyR 2 B()Q/c = — 2QP/c, so since near the star we need j3 to 
carry the torque away from the star we set /3 = —2Q(t)P/c 
for 2 « 0, P < P s and /3 = -2Q(t - z/c)P/c for z large 
since that gives the field found in Section 4. We therefore 
choose as our approximation /3 = — 2Q.(t — z/c)P/c whenever 
z > and P < Ps. For P > Ps, /3 = 0, see Fig 3. Since we 
have used the P of the non-rotating configuration which will 
in practice be modified by the rotating charges this is only 
approximate; however we note that in our rotating cylinder 
the B z remained at Bo even in the presence of the rotation 
so at large z the proposed solution is exact. 

Now we consider the case in which the light cylinder 
is of smaller radius than the asymptotic radius, a, of the 
separatrix of the static field. Close to the light cylinder the 
interstellar magnetic field is far too wea k to have any effect 
so for the aligned pulsar, the work of iGoldreich fe Julian! 




Figure 2. The magnetic field of a dipole embedded inside a uni- 
form external magnetic field of the same polarity. The thick red 
line is the separatrix of the field lines that are attached to the 
dipole field and those detached. Note that amongst the attached 
field lines the ones emerging near the poles reach infinity whereas 
the ones emerging from the area near the equator close in the 
other hemisphere of the dipole. 



( 1969) as modified by Scharlemann fc Wagoner] (|l973t ) and 
most recently treated iGoodwin et all (|2004T l is appropriate. 
However in those treatments the far field has a weak poloidal 
component and spirals off to infinity with that component 
not far from radial. There the field will be radically modified 
if the pulsar resides in an asymptotic field from the interstel- 
lar medium because asymptotically the field from the pulsar 
must align along the cylinders and produce Alfven waves 
propagating along them. QR will be much greater than c. 
Fig. 4 illustrates the type of the field configuration that we 
envisage. All the pulsar field lines that do not close within 
the light cylinder end up winding around the interstellar 
field and their Q is still that of the pulsar! 

The flux function for the pure dipolar field is P — 
2-kM{1 — fJ. 2 )/r so the flux that fails to reach the equator 
within the light cylinder in such a field would be P c = 
2ttMQ/c. In reality the field is distorted by the rotation to 
be more like P = 27rM/r(l + (l/2)(r/r c ) 3 )(l - /?) for r < 
r c = c/fi which gives P c = 3nMCl/c. A very rough model 
that gets the fluxes right (but is not force-free) is given by 
the P above for r < r c , P = P c (l — fi 2 ), for r c < r < r s where 
r 2 s = P c /(3nB ) and P = 7r£ r- 2 (l + 2(r 3 /r) 3 ) for r > r s . 
As before the toroidal field is given as in equation (28) with 
P(P) = -2fi(t - \z\/c)P/c for P < P c and zero elsewhere. 

To see the orders of magnitude consider a pulsar whose 
polar field is B p and suppose the radius of the polar cap 
whose field does not close within the light cylinder is R p and 
the inter-stellar field is Bo- Then the asymptotic radius of 
the field attached to the pulsar will be a = R p yj B p /Bo suit- 
able values might be B p = 3 x 10 12 G, B = 3 x 10" 6 G and 
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Figure 3. The magnetic field of a rotating dipole embedded in- 
side a uniform external magnetic field of the same polarity. The 
discontinuity of the toroidal component of the magnetic field on 
the separatrix leads to the formation of a current sheet. We have 
assumed that the combination of the pulsar magnetic field, in- 
sterstellar magnetic field and rotation velocity is such that the 
light cylinder occurs outside the co-rotating loops. It is possible 
that the light cylinder while being outside the co rotating loop 
intersects with the separatrix, in which case some of the pulsar's 
twisted field lines will pass through the light cylinder. 




Figure 4. The magnetic field of a rotating dipole embedded in- 
side a uniform external magnetic field of the same polarity. The 
discontinuity of the toroidal component of the magnetic field on 
the separatrix leads to the formation of a current sheet. Unlike 
Fig. 3 the combination of the fields and rotation velocity is such 
that the light cylinder occurs near the star. In this case a sheet 
current forms on the equator. 




Figure 5. The magnetic field of a dipole embedded inside a uni- 
form external magnetic field of the opposite polarity. The at- 
tached field lines to the dipole are enclosed within a spherical 
surface. The external field lines arc deformed from being uniform 
close to the equator. 




Figure 6. The magnetic field of a rotating dipole embedded in- 
side a uniform magnetic field of the opposite polarity. The field's 
configuration is similar to that of Fig. 4, however the direction of 
the field lines is different. There is a discontinuity on the sepa- 
ratrix of the poloidal component of the magnetic field enhancing 
the current sheet that already exists because of the discontinuity 
on the toroidal component of the field. 



R p — 1km giving a = 10 1 



the pulsars rotation rate Q and radius r s by R p — r s <J r s Q/c 
which gives R p — 1km for r s — 11km for a pulsar with a 
period of 28ms. More generally a = r s y (flr s /c)B p /Bo oc 
T _1//2 where T is the pulsar period. For the numbers above 
but general T, we find Ha = 3.5 x 10 3 c(T/s) _3/2 which is 
many times the speed of light. If the Alfven wave advances at 
c then in the interstellar gas B^/Bq = 3.5 x 10 3 (T/s)~ 3/2 in 
the idealised model and this is increased by the factor c/ca if 
the Alfven speed is ca - However such highly wound magnetic 
fields are likely to be highly unstable. In the ideal case the 
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plasma has to move along the field so that v z — (Bo/B^ca- 
The field winds through the plasma like an Archimidean 
screw winds up the water. While these estimates will be 
complicated by instabilities they do indicate that pulsars 
generate much enhanced fields in the interstellar medium 
albeit over tubes of force with diameters similar to that of 
the solar system. 

A different but interesting situation occurs when the 
interstellar field is aligned oppositely to magnetic moment 
of the pulsar. The total flux of the pulsar outside its light 
cylinder, P c , will be unable to counteract the flux of the 
interstellar field equatorially all the way to infinity since the 
latter grows like R 2 . Furthermore the pulsar's flux can not 
cross its equatorial plane outside the light cylinder as by 
symmetry it would have to cross perpendicularly and the 
plasma attached would then move faster than light! However 
the flux from the pulsar has to reach infinity since it carries 
the flow of plasma pushed by the Archimedian screw of the 
whirling field. The only place where the flux from the pulsar 
can deform the field over large distances is close to the axis, 
indeed in the non-rotating case there would be a neutral 
point on axis (jDungeyl 19581 ). see Fig. 5. As the pulsar rotates 
the equatorial diameter of the field structure around it grows 
and this neutral point is sent to infinity along the axis. An 
oppositely directed poloidal field is generated near the axis 
and is accompanied by the B$ and the plasma flow due 
to the pulsar's rotation. We illustrate both the weak spin 
case when all the star's field lies within the light cylinder 
and rotates with the star Fig. 5 and the pulsar case just 
described in Fig. 6. Corresponding to any solution such as 
that shown in Fig. 4. there is an other solution with the 
field reversed outside the separatrix. If the longitudinal field 
inside the separatrix is small compared to the toroidal field 
there, the field pressures and tension still balance along the 
separatrix. The argument given above leads us to believe 
that this is the field configuration of the steady anti-aligned 
pulsar. 



to the fire-hose instability. However there remains the in- 
teresting possibility that the magnetic fields generated over 
small regions by the thrashing field lines might be surpris- 
ingly large although smaller than the maxima estimated via 
these axially symmetric models. 
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6 CONCLUSIONS 

In this paper we have studied analytical solutions of rotating 
magnetic fields emerging from dipoles or discs embedded in 
uniform magnetic fields. This study is based on the assump- 
tion that the dipole field is either aligned or counter-aligned 
to the external field. Fields at an angle i to the dipole may 
be considered by taking the solutions above with So re- 
placed by Bo cos i and superposing a lateral field parallel to 
the plane of the disk of magnitude Bosini, but this gives 
gives only orders of magnitude and is no substitute for a full 
non-axially symmetric treatment. 

We find that the spinning field will form Alfven waves in 
the interstellar gas. There is no need to confine these Alfven 
waves inside the light cylinder, as the plasma can have a 
velocity component parallel to the magnetic field and not 
exceed the speed of light and we find that for typical pulsar 
values they will occupy cylinders whose cross sectional areas 
are of the size of the solar system. These field configurations 
with tightly wound fields outside the light cylinder are prob- 
ably very unstable, in particular the pressure of the whole 
tube along the axis is very large so any departure from ax- 
ial symmetry will be enhanced by an instability analogous 



